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Pairing in high-density neutron matter including short- and long-range correlations 
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The influence of short-range correlations (SRC) on the spectral distribution of neutrons is in¬ 
corporated in the solution of the gap equation for the 3 P 2 — 3 F 2 coupled channel in pure neutron 
matter at high density. This effect is studied for three different realistic interactions. The gap in 
this channel is strongly suppressed by these correlations but does not vanish. For a consistent treat¬ 
ment we also include for the first time the effect of long-range correlations (LRC) by incorporating 
polarization terms in addition to the bare interaction. This allows the neutrons to exchange density 
and spin fluctuations governed by the strength of Landau parameters with values that are consistent 
with the available literature. While these LRC have an antiscreening tendency, they only slightly 
increase the gap in the 3 P 2 — 3 F 2 coupled channel for all three realistic interactions as long as SRC 
are included. All three interactions generate maximum gaps around 0.1 to 0.2 MeV at most with 
a small dependence on the hardness of the interaction. These results are relevant for the cooling 
scenarios of neutron stars, in particular the young neutron star in Cassiopeia A. 

PACS numbers: 21.65.-f, 21.65.Cd, 24.10.Cn, 26.60.-c, 74.20.Fg 


Pulsar glitches are considered prime evidence of super¬ 
fluid properties in neutron stars on account of their long 
relaxation times Eli- The cooling scenario of a neutron 
star also depends sensitively on the possible superfluid¬ 
ity of neutrons and superconductivity of protons in the 
interior of the system |3[ . Recent observations of a young 
neutron star in Cassiopeia A Q and its rapid cooling [5j 
have rekindled interest in high-density neutron superflu¬ 
idity @-jH- Several scenarios with apparently different 
mechanisms appear to be able to explain the observed 
cooling observation @-H] but all require a contribution 
of neutron pairing with 3 P 2 — 3 F 2 quantum numbers. 
Longer observation of this system may shed further light 
on this issue. Meanwhile, it appears worthwhile to con¬ 
sider the present status of high-density neutron pairing. 


At low density there is universal agreement that neu¬ 
tron pairs with 1 S'o quantum numbers will be super¬ 
fluid on account of the exhibited attraction of the 

corresponding low-energy phase shift. At higher energy, 
the P -wave phase shift of the coupled 3 P 2 — 3 F 2 channel 
exhibits more attraction and therefore becomes a candi¬ 
date for neutron pairing at higher density While 1 5'o 
gaps in standard BCS calculations are quite similar for 
different realistic interactions, there is more variation in 
the gap obtained for the 3 P 2 — 3 F 2 channel. Depending 
on the choice of the nucleon-nucleon (NN) interaction, 
the results range from about 0.5 to 0.9 MeV for the max¬ 
imum gap with a density range associated with Fermi 

momenta from 1.8 to 2.5 fm -1 a eh eh. 


Another example of an attractive interaction pertain¬ 
ing to a coupled channel occurs in symmetric nuclear 
matter for isoscalar pairs with 3 Si- 3 Z?i quantum num¬ 


bers. We build on the experience with pairing in this 
coupled channel to motivate our treatment of neutron 
3 P 2 — 3 F 2 pairing. We note that early calculations 
around normal nuclear matter density generated a sizable 
gap of around 10 MeV (Tfll [771 with 3 Si~ 3 Di quantum 
numbers when employing e.g. Brueckner-Hartree-Fock 
(BHF) single-particle (sp) energies. The empirical data 
of finite nuclei give no indication for such strong proton- 
neutron pairing correlations with corresponding gaps as 
large as 10 MeV. It is therefore plausible that the evalu¬ 
ation of the 3 Si~ 3 Di gaps in infinite matter according to 
Refs. [l3l - fl7| suffers from the lack of inclusion of relevant 
physics that can suppress such a large gap. 

Nucleons in nuclei are well documented to behave dif¬ 
ferently than expected on the basis of a pure indepen¬ 
dent particle model (IPM). Additional correlations be¬ 
yond the Pauli principle of the IPM are clearly estab¬ 
lished in the analysis of the (e, e'p) reaction of closed-shell 
nuclei EH leading to removal probabilities for valence 
protons that are reduced by about 0.35 from the IPM 
value of 1. A substantial contribution to this reduction, 
about 0.10 to 0.15, is a consequence of the depletion of 
the Fermi sea due to short-range correlations (SRC) EH- 
The complementary admixture of a modest amount of 
high-momentum components in the nuclear ground state 
has been experimentally observed (20j. 

Such correlations can be accurately and completely 
treated in infinite nuclear systems where the conse¬ 
quences of summing ladder diagrams can be included into 
the self-energy and off-shell propagation can be incorpo¬ 
rated self-consistently EH- While such a procedure still 
encounters some numerical issues at zero temperature no 
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such problems are encountered at sufficiently high tem¬ 
perature 21 - 24]. At finite temperature it is thus possible 
to fully account for the influence of SRC on the propa¬ 
gation properties of nucleons in symmetric nuclear and 
neutron matter. 

At T — 0 ladder-diagram summations can lead to pair¬ 
ing instabilities when the possibility of anomalous (pair) 
propagation is not included. Since the influence of pair¬ 
ing on the normal nucleon propagator is confined to en¬ 
ergies around the Fermi energy that are of the order of 
the size of the gap, one may reasonably assume that the 
normal self-energy will hardly be affected by pairing cor¬ 
relations J2fj. Making this assumption and extrapolat¬ 
ing the temperature dependence of the self-consistently 
calculated normal self-energy to lower temperatures in¬ 
cluding T = 0, it was shown in Ref. [2(| that due to the 
influence of SRC the 3 Si~ 3 D 1 gap at normal density van¬ 
ishes, reconciling empirical information with many-body 
calculations of pairing in nuclear matter. We therefore 
follow such a procedure for the calculation of the influ¬ 
ence of SRC on the pairing properties in neutron matter 
at high density for the 3 P 2 — 3 F 2 coupled channel. 

Another important ingredient for the calculation of 
pairing correlations at high density in neutron matter is 
the modification of the pairing interaction when nucleons 
exchange the low-energy, possibly collective, excitations 
of the medium. At small momentum transfer such excita¬ 
tions are determined by the Landau parameters that gov¬ 
ern the density and spin excitations of the medium [241. 
The first correction to be added to the bare NN in¬ 
teraction is given by the exchange of one particle-hole 
bubble between the neutrons. This mechanism can be 
extended to exchanging an infinite-order random-phase- 
approximation (RPA) bubble series. Such contributions 
arise since the gap equation itself generates the contribu¬ 
tion of ladder diagrams [3, [H, S3] ■ The accurate calcula¬ 
tion of Landau parameters depends on many ingredients 
some of which were discussed in Ref. 0]. We will also 
include the contribution of these so-called polarization 
terms that reflect LRC and follow the procedure outlined 
in Ref. [2f| for neutron 1 So superfluidity and extend it to 
the case of the 3 P 2 — 3 F 2 coupled partial waves thereby 
generating for the first time a complete treatment of both 
SRC and LRC for this pairing channel. 

We start with a brief and schematic reminder of the 
self-consistent summation of ladder diagrams that was 
used for the construction of the ingredients of the pair¬ 
ing calculation jdO;, 311. The proper treatment of SRC 
and tensor correlations in the nuclear medium is ac¬ 
complished by constructing the in-medium Lippman- 
Schwinger equation, schematically represented by: 


r = v 


vg° n T. 


(i) 


The NN interaction is denoted by V and the noninteract¬ 
ing but dressed two-particle propagator is given by the 


convolution of two dressed sp propagators 

g° u {k,k'-,n) = J ^-^-A{k,u)A{k',u') 


1 - f{u) - /(£*/) 


— cu — ( 
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where f(u>) = (1 + _1 is the Fermi-Dirac distri¬ 

bution and A the sp spectral function. The temperature, 
T, is applied as an external condition. The chemical po¬ 
tential, fx, is obtained from the normalization to the den¬ 
sity of the momentum distribution. Self-consistency is 
imposed all the way through in our calculations. The 
sp spectral functions that enter Eq. ©, for instance, 
are obtained from the T— matrix itself. In the ladder 
approximation, this effective interaction determines the 
imaginary part of the self-energy Si- 

The dispersive contribution to the real part of the self¬ 
energy can be obtained from a dispersion relation 24]. In 
addition the energy-independent correlated Hartree-Fock 
contribution involving the actual momentum distribution 
must be included as well. From Dyson’s equation we then 
obtain the sp spectral function: 


A(k,w) = 


— 2ImE(/c, w) 


[w - & - Re£(fc,cc)]2 + pmE(fc,w)P ’ 


(3) 


which closes the self-consistency loop that treats all parti¬ 
cles in the same manner, providing appropriate feedback 
to the different ingredients of the calculations. 

The proper treatment of pairing requires the introduc- 

25] with correspond- 


tion of anomalous propagators [24], 
ing anomalous self-energy terms in so-called Gorkov 
equations ]25j, 32, S3] that couple normal and abnor¬ 
mal (pair) propagators. The resulting expression for the 
anomalous self-energy A can be written in the form of the 
usual gap equation in a partial wave expansion [26l. l34| 


,JST 




(4) 


For a coupled channel we employ the spherical gap 


,JST( 


(*) = [(A/fJ.^fc)) 2 + (A /JJ +1 (*)) ! 


1/2 


(5) 


in E k = + ]A JST (fc)] 2 . With Xfc = £t - £fc, while e k 

represents the kinetic energy or medium-modified spec¬ 
trum without pairing, one obtains the conventional BCS 
gap. The inclusion of SRC requires the calculation of 
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where A s is the spectral function that includes a nontriv¬ 
ial pairing solution to the gap equation. It is therefore 
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appropriate to consider Eq. © as a generalization of the 
usual gap equation when Ek is obtained from Eq. © , as 
it accounts for the spreading of sp strength due to SRC 
and tensor correlations !H • For practical calcula¬ 
tions we assume that different partial-wave projections 
do not couple through the angular dependence of the in¬ 
teraction and the propagators. We therefore employ the 
spherical version of the quasiparticle energy as in Eq. ©• 

The essential calculation that provides the input for 
the solution of the gap equation of Eq. © is provided by 
a proper temperature extrapolation [361 ] of the convolu¬ 
tion of spectral functions generated by the self-consistent 
normal self-energy in which A s in Eq. © is replaced by 
A whereby we obtain a generalized version \k instead of 
Xk [HI at a given temperature. For the present work we 
have explored an expansion of both the real and imag¬ 
inary part of the normal self-energy in even powers of 
the temperature. This expansion is constrained by rele¬ 
vant macroscopic thermodynamical quantities. The zero- 
temperature limit can then be used to construct the rele¬ 
vant spectral functions modified by SRC. This procedure 
was tested by comparing with the results of Ref. [26] for 
nuclear and neutron matter in which a slightly different 
extrapolation method was employed. 

We also report results for the inclusion of polarization 
terms in the pairing interaction. We have followed the 
procedure proposed in Ref. [2!| and generalized it to the 
case of the 3 P 2 — 3 F 2 coupled channel. The physical ex¬ 
change of low-lying possibly collective density and spin 
modes can be included by adding the corresponding in¬ 
teraction terms to the bare NN interaction. These excita¬ 
tions are dominated by the strength of the corresponding 
Landau parameters which govern the low-lying spectrum 
of the fluid. The possibility that the presence of the pion- 
exchange tensor interaction strongly influences the spin 
mode [371 ] cannot be excluded but is beyond the scope 
of the present work jslj]. The corresponding interaction 
that treats LRC for the 1 S'o channel is given by 


is attractive and the second term is repulsive but domi¬ 
nates on account of the spin factor ultimately leading to 
additional suppression of the gap in this channel. 

For the 3 P 2 — 3 F 2 channel which involves spin-1 pairs, 
the sampling over density and spin modes yields 


1 


1 


Vlrc (?) — ^Gp h Gp h A n (q) + -Gp h Gp h A 1 (q), (9) 

with both terms yielding attraction. Contrary to the 1 iS'q 
channel this contribution leads to antiscreening of the gap 
as it represents an attractive interaction. The Landau 
parameters C s from Ref. [29| are extrapolated to higher 
densities in a smooth way to approximately fulfill the 
forward scattering sum rule for these parameters [HI 
leading to their sign reversal at a density corresponding 
to about 1.8 fm _1 . 

Results for the 1 5' 0 channel including SRC are con¬ 
sistent with the results of Ref. [26} for the charge- 
dependent Bonn (CDBonn) interaction [40] and Argonne 
V18 (AV18) interaction of Ref. 411. The inclusion of po¬ 
larization terms on top of the bare interaction also re¬ 
produce the results of Ref. [HI ■ The main effect of the 
inclusion of SRC in this channel is to reduce the value 
of the gap by about 0.75 MeV compared to the BCS 
results that employ a free sp spectrum and generate a 
gap closure at a corresponding lower density. We note 
that the inclusion of self-energy effects due to SRC re¬ 
moves sp strength from the vicinity of the Fermi energy 
and distributes it far and wide without compensating the 
loss of strength near the Fermi energy. This feature has 
the inescapable consequence that a pairing solution is 
less likely and therefore the gap is suitably reduced or 
vanishes. Other calculations that take such effects into 
account yield identical conclusions J42[j4hj]. 

Results for the Idaho next-to-next-to-next-to-leading- 
order (N3LO) chiral-perturbation-theory potential of 
Ref. |46[ have also been generated. Our results for the 
three realistic interactions including both SRC and LRC 


Vlrc(i) — -^G°p h G 0 ph A a (q) - -Gp h Gp h A 1 (q), (7) 

C< 

where G p J[ h are particle-hole transformed G-matrix ele¬ 
ments averaged around the Fermi energy [29| . The iter¬ 
ated bubble series is represented by 


A S (q) 


A 0 (g) 

1 - A 0 (q) * C s ' 


( 8 ) 


The C s correspond to Landau parameters. The density 
mode with total spin 0 in the particle-hole channel is de¬ 
termined by C° which is attractive at low density and 
the spin mode with total spin 1 by C 1 which is repulsive 
but has similar magnitude, while both exhibit a modest 
density dependence [29j. The static Lindhard function 
Ao (q) is iterated to all orders in this procedure to gener¬ 
ate the A s which are negative. The first term in Eq. © 



FIG. 1. Results for the 1 S'o gap in neutron matter for three 
different realistic interaction with the inclusion of SRC and 
separately with additional polarization terms. 
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FIG. 2. Gap for 3 P 2 — 3 P 2 neutron pairing as a function of the Fermi momentum in pure neutron matter using the AV18 (left 
panel), CDBonn (middle panel), and N3LO interaction (right panel). For symbols see the text. 


are very similar and characterized by a maximum gap 
of about 1.8 MeV at a density corresponding to a Fermi 
momentum of a little more than 0.8 fm _1 . These results 
appear very robust and confirm the general understand¬ 
ing of the gap in this channel in terms of its size and 
density dependence 0 - Results are shown in Fig. [T] for 
the CDBonn interaction by the dash-dot (SRC only) and 
long-dashed curve (SRC + polarization), N3LO curves 
by the short-dashed (SRC) and solid curve (SRC + po¬ 
larization), and AV18 curves by the dash-dot-dot (SRC) 
and dash-dah-dot curve (SRC + polarization). 

In the case of the 3 P 2 ~ 3 F 2 channel, the spin recoupling 
from particle-hole to particle-particle leads to an inter¬ 
action that enhances the pairing gap and can therefore 
be characterized as antiscreenin g; a ccording to Eq. ©■ 
Preliminary results reported in [471 ] generated no pair¬ 
ing solution in this channel. A more detailed treatment 
of the extrapolation to T = 0 has recently been imple¬ 
mented [0 yielding a small but nonvanishing gap for all 
three interactions. We show in Fig.[2]the results for these 
interactions including the standard BCS result (solid 
line), additional inclusion of polarization (dashed), in¬ 
clusion of SRC (open circles), and inclusion of both SRC 
and LRC (full circles). A fit to the calculated data points 
is represented by solid curves using the parametrization 
of Ref. 3 (see the table in Ref. [36|). The calculations 
for the N3LO interaction have been restricted to densities 
somewhat below its 500 MeV cut-off 481. 


The BCS results for these interactions are quite simi¬ 
lar up to a density corresponding to fcp « 2.0 fm -1 but 
slightly differ in their maximum gap and the density for 
gap closure. CDBonn generates the largest values for 
these quantities. The effect of adding only polarization 
to the BCS treatment then yields dramatically different 
results illustrated by the dashed curves in Fig. [2] with 
maxima (not shown) corresponding to 1.46 (AV18) , 2.68 
(CDBonn), and 1.33 MeV (N3LO), respectively. When 
only the SRC are incorporated the results for the three 
interactions indicated by open circles become much closer 


with maximum gaps ranging from 0.04 MeV (AV18), 0.07 
(CDBonn) to a about 0.1 MeV (N3LO). The size of the 
employed symbols in the figures reflect a reasonable es¬ 
timate of the error associated with the T = 0 extrap¬ 
olation filil ]. The additional inclusion of LRC to the 
SRC effect enhances the gap slightly for all three interac¬ 
tions with the softest interaction (N3LO) generating the 
largest gap of about 0.2 MeV, CDBonn 0.1 MeV, and the 
harder AV18 a maximum of about 0.08 MeV. 

We conclude that the inclusion of SRC has a dramatic 
but not unexpected effect on the solution of the gap equa¬ 
tion in high-density neutron matter for the 3 P 2 — 3 F 2 cou¬ 
pled channel. This conclusion is valid for the hard AV18, 
the softer CDBonn, as well as the very soft N3LO inter¬ 
action. Our results indicate that a small gap remains at 
high density irrespective of the realistic interaction. The 
treatment of LRC suggests that SRC generate the dom¬ 
inant effect but they do enhance the gap to maximum 
values ranging from 0.1 to 0.2 MeV at densities corre¬ 
sponding to Icf ~ 1-7 fm -1 and kF ~ 1-9 fm -1 , respec¬ 


tively. More sophisticated treatments of LRC m appear 
unlikely to change these results by much as the effect of 
SRC remains dominant. We note that the smallness of 
the gaps may require an accuracy that may be difficult 
to reach by Monte Carlo calculations in the near future. 
The inclusion of three-body forces associated with the 
N3LO interaction appears possible 50|, 51] but may not 
change these conclusions according to Ref. [521. 

The present results yield for the first time a consistent 
treatment of both SRC and LRC for the 3 P 2 — 3 F 2 cou¬ 
pled channel and suggest that a small gap in a narrow 
density regime can be expected at densities where real¬ 
istic interactions are still constrained by NN scattering 
data. The impact of our conclusions on cooling calcula¬ 
tions is likely to be qualitatively significant. A complete 
understanding of cooling phenomena requires a consis¬ 
tent description of both the equation of state as well as 
pairing properties. The present effort represents a first 
step in providing these ingredients. 


















5 


This material is based upon work supported by 
the U.S. National Science Foundation under grant 
PHY-1304242, Consolider Ingenio 2010 Programme 
CPAN CSD2007-00042, Grant No. FIS2011-24154 from 
MICINN (Spain), and Grant No. 2014SGR-401 from 
Generalitat de Catalunya (Spain); by STFC, through 
GrantsNo. ST/1005528/1 and No. ST/J000051; and by 
NewCompstar, COST Action MP1304. 


[1] A. B. Migdal, Nucl. Phys. 13 , 655 (1959). 

[2] N. Charnel and P. Haensel, Living Rev. Relativity 11 , 10 
(2008). 

[3] D. Page, J. M. Lattimer, M. Prakash, and A. W. Steiner, 
Astrophys. J. Suppl. Ser. 155, 623 (2004). 

[4] W. C. G. Ho and C. O. Heinke, Nature 492 , 71 (2009). 

[5] P. S. Shternin, D. G. Yakovlev, C. O. Heinke, W. C. G. 
Ho, and D. J. Patnaude, Mom. Not. R. Astron. Soc. 
412 , L108 (2011). 

[6] D. Page, M. Prakash, J. M. Lattimer, and A. W. Steiner, 
Phys. Rev. Lett. 106 , 081101 (2011). 

[7] D. Blaschke, H. Grigorian, and D. N. Voskresensky, 
Phys. Rev. C 88, 065805 (2013). 

[8] W. C. G. Ho, K. G. Elshamouty, C. O. Heinke, and A. Y. 
Potekhin, Phys. Rev. C 91 , 015806 (2015). 

[9] D. J. Dean and M. Hjorth-Jensen, Rev. Mod. Phys. 75, 
607 (2003). 

[10] S. Gandolfi, A. Y. Illarionov, F. Pederiva, K. E. Schmidt, 
and S. Fantoni, Phys. Rev. C 80 , 045802 (2009). 

[11] M. Baldo, O. Elgaroy, L. Engvik, M. Hjorth-Jensen, and 
H. J. Schulze, Phys. Rev. C 58 , 1921 (1998). 

[12] V. A. Khodel, J. W. Clark, M. Takano, and M. V. 
Zverev, Phys. Rev. Lett. 93, 151101 (2004). 

[13] T. Aim, G. Ropke, and M. Schmidt, Z. Phys. A 337, 
355 (1990). 

[14] B. E. Vonderfecht, C. C. Gearhart, W. H. Dickhoff, 
A. Polls, and A. Ramos, Phys. Lett. B 253, 1 (1991). 

[151 M. Baldo, I. Bombaci, and U. Lombardo, Phys. Lett. B 
283 , 8 (1992). 

[16] T. Takatsuka and R. Tamagaki, Suppl. Prog. Theor. 
Phys. 112, 27 (1993). 

[17] M. Baldo, U. Lombardo, and P. Schuc.k, Phys. Rev. C 
52 , 975 (1995). 

[18] L. Lapikas, Nucl. Phys. A553, 297c (1993) 

[19] W. H. Dickhoff and C. Barbieri, 
Prog. Part. Nucl. Phys. 52 , 377 (2004) 

[20] D. Rohe et al, Phys. Rev. Lett. 93, 182501 (2004) 

[21] E. P. Roth, Self-consistent Green’s Functions in Nuclear 
Matter, Ph.D. thesis, Washington University in St. Louis 
( 2000 ). 

[22] T. Frick, Self-consistent Green’s Functions in Nuclear 
Matter at Finite Temperature , Ph.D. thesis, University 


of Tubingen (2004). 

[23] A. Rios, Thermodynamical Properties of Nuclear Matter 
from a Self-Consistent Green’s Function Approach, Ph.D. 
thesis, University of Barcelona (2007). 

[24] W. H. Dickhoff and D. Van Neck, Many-body theory ex¬ 
posed!, 2nd ed. (World Scientific, 2008). 

[25] A. B. Migdal, Theory of Finite Fermi Systems (Inter¬ 
science, New York, 1967). 

[26] H. Miither and W. H. Dickhoff, Phys. Rev. C 72, 054313 
(2005). 

[27] W. H. Dickhoff, Phys. Lett. B 210, 15 (1988). 

[28] W. H. Dickhoff and H. Miither, Nucl. Phys. A 473, 394 
(1987). 

[29] L. G. Cao, U. Lombardo, and P. Schuck, Phys. Rev. C 
74, 064301 (2006). 

[30] A. Rios, A. Polls, and W. H. Dickhoff, 
Phys. Rev. C 79, 064308 (2009) 

]31] A. Rios, A. Polls, and W. H. Dickhoff, Phys. Rev. C 89, 
044303 (2014). 

[32] L. P. Gorkov, Sov. Phys.-JETP 7, 505 (1958). 

[33] V. Soma, T. Duguet, and C. Barbieri, Phys. Rev. C 84, 
064317 (2011). 

[34] P. Bozek, Nucl. Phys. A 657, 187 (1999). 

[35] Y. Dewulf, W. H. Dickhoff, D. Van Neck, E. R. Stoddard, 
and M. Waroquier, Phys. Rev. Lett. 90, 152501 (2003) 

[36] See Supplemental Material at 

http://link.aps.org/supplemental/10.1103/PhysRevLett. 

[37] W. H. Dickhoff, A. Faessler, and H. Miither, Phys. Rev. 
Lett. 49, 1902 (1982). 

[38] S. S. Pankratov, M. Baldo, and E. E. Saperstein, Phys. 
Rev. C 91, 015802 (2015). 

[39] B. L. Friman and A. K. Dhar, Phys. Lett. B 85, 1 (1979). 

[40] R. Machleidt, F. Sammarruca, and Y. Song, 
Phys. Rev. C 53, R1483 (1996) 

[41] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, 
Phys. Rev. C 51, 38 (1995) 

[42] P. Bozek, Phys. Rev. C 62, 054316 (2000). 

[43] M. Baldo and A. Grasso, Phys. Lett. B 485, 115 (2000). 

[44] M. Baldo, U. Lombardo, H. J. Schulze, and Z. Wei, Phys. 
Rev. C 66, 054304 (2002). 

[45] P. Bozek, Phys. Lett. B 551, 93 (2003). 

[46] D. R. Entem and R. Machleidt, 
Phys. Rev. C 68, 041001 (2003) 

[47] D. Ding, S. J. Witte, W. H. Dickhoff, H. Dussan, A. Rios, 
and A. Polls, AIP Conf. Proc. 1619, 73 (2014). 

[48] S. Maurizio, J. W. Holt, and P. Finelli, Phys. Rev. C 
90, 044003 (2014). 

[49] A. Sedrakian, Phys. Rev. C 68, 065805 (2003). 

]50] A. Carbone, A. Polls, and A. Rios, Phys. Rev. C 88, 
044302 (2013). 

[51] A. Carbone, Self-consistent Green’s Functions with 
Three-body Forces, Ph.D. thesis, University of Barcelona 
(2014). 

[52] J. M. Dong, U. Lombardo, and W. Zuo, Phys. Rev. C 
87, 062801 (R) (2013). 





